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ABSTRACT 

We demonstrate that the high isotropy of the Cosmic Microwave Background 
(CMB), combined with the Copernican principle, is not sufficient to prove homogeneity 
of the universe - in contrast to previous results on this subject. The crucial additional 
factor not included in earlier work is the acceleration of the fundamental observers. 
We find the complete class of irrotational perfect fluid spacetimes admitting an 
exactly isotropic radiation field for every fundamental observer and show that they 
are Friedman-Lemaitre-Robertson- Walker (FLRW) if and only if the acceleration is 
zero. While inhomogeneous in general, these spacetimes all possess three-dimensional 
symmetry groups, from which it follows that they also admit a thermodynamic 
interpretation. 

In addition to perfect fluid models we also consider multi-component fluids containing 
non-interacting radiation, dust and a quintessential scalar field or cosmological constant 
in which the radiation is isotropic for the geodesic (dust) observers. It is shown that the 
non-acceleration of the fundamental observers forces these spacetimes to be FLRW. 

While it is plausible that fundamental observers (galaxies) in the real universe 
follow geodesies, it is strictly necessary to determine this from local observations for the 
cosmological principle to be more than an assumption. We discuss how observations 
may be used to test this. 



1. Introduction 

The high isotropy of the CMB is usuaUy taken as strong evidence that the universe is 
homogeneous and isotropic, i.e., is well described by an FLRW model. The principle justification for 
this is an important theorem of Ehlers, Geren and Sachs (1968) (based on earlier work by Tauber 
and Weinberg 1961), which states that if all observers in an expanding, dust universe measure an 
isotropic CMB then the universe is FLRW and the cosmological principle is valid. The importance 
of this theorem lies in the fact that it permits the homogeneity and isotropy of the universe to be 
deduced not from measurements of the actual isotropy of the universe about us, but from only 
measurements of the CMB, combined with the Copernican principle (that is, the assumption that 
all observers in the universe see the same degree of isotropy) . The Copernican principle is often 
regarded as a powerful but untestable assumption in cosmology, although there arc suggestions 
that it may be testable using the Sunyaev-Zeldovich effect, for example (Goodman 1995). Here we 
simply assume that the Copernican principle is valid and study the consequences of applying it to 
the observed high degree of isotropy of the CMB. That is, we examine spacetimes with an isotropic 
CMB for all observers. 

The EGS theorem has been generalised by Treciokas and Ellis (1971) to include an isotropic 
collision term. Ferrando, Morales, and Portilla (1992) find the general form of the energy-momentum 
tensor and Einstein's equations for spacetimes with an isotropic radiation field, and consider 
some special cases with anisotropic pressure. It has also been shown by Stoeger, Maartens and 
Ellis (1995) that the EGS theorem almost holds when applied to an almost isotropic radiation field. 

There are counterexamples to the spirit of the EGS theorem (that is, when some of the 
assumptions are relaxed the result fails to hold). In particular, Ellis, Maartens and Nel (1978) show 
that the result does not hold if the expansion is zero (which is obviously not relevant to cosmology) , 
and Ferrando et al. (1992) emphasise that homogeneity does not follow if there is anisotropic 
pressure in the energy-momentum tensor. Most importantly for the work presented here, though, 
is the result of Barrett and Clarkson (1999a), which shows that when the assumption of geodesic 
observers is relaxed there exist inhomogcneous perfect fluid (or scalar field) cosmologies with an 
isotropic CMB. (In fact, we show in this paper that the Stephani models considered in Barrett and 
Clarkson (1999a) are representatives of the only perfect fluid spacetimes that admit an isotropic 
CMB for all fundamental observers.) Nilsson et al. (1999) provide a counterexample to the almost 
EGS result when the Weyl curvature is not negligible. 

The basis of the EGS theorem is the Liouville equation for photons, which tells us that if a 
radiation field (i.e., a solution of the Liouville equation) exists such that for every observer on some 
timelikc congruence the radiation field is isotropic, then that congruence is (parallel to) a conformal 
Killing vector (CKV). This may be expressed more formally as follows (Ehlers et al. 1968; Ferrando 
et al. 1992): 

Theorem 1 A spacetime will admit an isotropic radiation field if and only if it is conformal to a 
stationary spacetime, which happens if and only if there is a velocity field w" satisfying 

<Jab = 0, (1) 

V[, (iib] - i^iifc]) - 0, (2) 

where dab, u°' and 6 are the shear, acceleration and expansion of u"" , respectively. 



(Then u° is the velocity field relative to which the radiation is isotropic, and is parallel to the 
CKV.) 

In the absence of some statement about the matter content of a spacetime, or further 
assumptions about the congruence u", Theorem n^ is all that can be said. In a cosmological 
context the simplest, and most common, assumption is that the matter is dust (implying that m" is 
geodesic), which leads to: 

Theorem 2 (Ehlers, Geren, and Sachs) // the fundamental observers in a dust spacetime see 
an isotropic radiation field, then the spacetime is locally FLRW. 

Alternatively, we can simply assume that u" is geodesic. The existence of an isotropic 
radiation field then ensures (for non-zero expansion) that the energy flux relative to m° is zero. 
If the anisotropic stress tensor is zero at any instant (so that the energy-momentum tensor has 
perfect fluid form) then it will remain zero and the spacetime will be FLRW (Ferrando et al. 1992, 
Corollary 1; but note that their statement that the anisotropic stress is invariant along u° in general 
is misleading - from Eqs. (31) and (40) of Ellis 1998 we have ■7r(^ab) c« 0-Kab)- 

It is worth emphasising that in applications of the above results to cosmology the motion of 
the fundamental observers must be identified with the congruence w". For example, in §2.2| all 
Stephani models are conformally flat, and therefore conformally stationary, but for most of these 
spacetimes the fluid congruence is not aligned with the timelike CKV. 

The matter content of the universe is not precisely known. Certainly, there is a large number 
of possible contributors, including hot and cold dark matter (in their various manifestations), 
electromagnetic fields etc., as well as the more obvious radiation and baryonic matter. In particular, 
the type la supernova results of Perlmutter et al. (1999) suggest that an important component 
may be a 'quintessential' scalar field. However, the forms of matter that are thought to contribute 
significantly to the energy-momentum tensor may be treated in general as perfect fluids. That is, 
their energy-momentum tensor may be written in the form 

Tab ^ IJ-UaUb + phab, (3) 

where /i and p are the energy density and pressure, u"^ is the timelike velocity congruence of 
the fluid, and hab is the spatial projection tensor associated with u". Scalar fields may also be 
written in this form, with u" parallel to the gradient of the scalar field, provided that this is 
timelike (see §0). Note that if several such components are present there is no reason why their 
fundamental congruences (the m°'s) should be parallel. If they are, then they may be treated as 
effectively a single perfect fluid with the energy densities and pressures added together. If not, the 
decomposition of the energy-momentum with respect to u° for one such fluid (as in equation (14) of 
Ellis 1998) will contain energy flux and anisotropic stress terms from the other fluids (again, see Sg). 
The fundamental observers will be associated with one such congruence. Usually the fundamental 
observers in standard models of the universe are associated with a dust-like (p — 0) component, 
with the result that the acceleration ii"' of the fundamental congruence is zero. However, we wish 
to study the consequences of relaxing this assumption and consider models with acceleration. 
This acceleration must be caused by some non-gravitational force (typically pressure gradients for 
perfect fluid spacetimes, but in principle it could be the result of a coupling between the fluid and 
some other component such as the electromagnetic field). 



With this in mind, in this paper we consider perhaps the two simplest generahsations of the 
dust hypothesis. Firstly we imagine that the dominant form of matter is a single irrotational 
perfect fluid. We do not specify what form of matter this corresponds to, but we allow pressure 
gradients that give rise to acceleration. Secondly we consider cosmological models in which more 
than one matter component makes a significant contribution to the energy density and dynamics 
of the universe. Specifically we consider 'QCDM' models, containing a non-interacting mixture of 
radiation, dust (CDM) and a scalar field (or a cosmological constant). The observers are associated 
with the CDM component, and are therefore geodesic. The difference between this and other 
theorems assuming geodesic observers is that the the scalar field component can introduce effective 
energy flux and anisotropic stresses relative to the dust congruence, and so the matter need not 
behave as a perfect fluid. 

In the following section we flnd all irrotational perfect fluid solutions admitting an isotropic 
radiation field for the fundamental observers, showing that they form a subclass of the Stephani 
spacetimes and are FLRW if and only if the acceleration vanishes. Then in §0 we examine QCDM 
models and prove that such models must be homogeneous and isotropic if they admit an isotropic 
radiation field. Finally, in sy we emphasise the importance of acceleration for these results and show 
that the acceleration of the fundamental congruence is, in principle, detectable, and measureable. 



in galaxy surveys. Two appendices contain results relating to §2.2 



2. The Irrotational Perfect Fluid Solutions. 

We wish to consider the constraints imposed by the existence of an isotropic radiation 
field for the fundamental observers on perfect fluid spacetimes in which the rotation of the 
fundamental congruence is zero. Since it follows from Theorem |l| that the shear of the fundamental 
congruence must also be zero we immediately know that all the acceptable solutions are members 
of the Stephani-Barnes family, which is the family of all shear-free, irrotational, expanding (or 
contracting) perfect fluids (see Barnes 1973; Krasinski 1989, 1997). It only remains, then, to impose 
condition (ph of Theorem H and thus find the sub-class of Stephani-Barnes models which admit an 
isotropic radiation field for all fundamental observers. 

The Stephani-Barnes family contains the Barnes solutions, which are of Petrov type D, and the 
Stephani models, which are conformally flat, although these two classes overlap where the Barnes 
solutions degenerate to type O (these solutions then become Stephani models with symmetry) . The 
FLRW models are a subcase of these solutions. The Barnes spacetimes all possess symmetry (see 
below), whereas the Stephani spacetimes, in general, do not. In all cases the metric in comoving 
coordinates can be written in the form (see Krasinski 1989; although here we use the same notation 
as Krasinski 1997): 

ds^ = V-"^ {-{FVtfdt^ + dx^ + di/ + dz^} (4) 

where F — F[t) and V — V{t^ x, y, z). F(t) is arbitrary, but there are some restrictions on the form 
of V depending on the symmetries of the solution, and these will be discussed in due course (but 
the impatient reader may wish to note Eqs. ^ O, O, and 11^). Expressions for the energy density 
and pressure can be found in Krasiiiski (1997). 



The fluid velocity is given by (without loss of generality we can assume that V > 0) 

V 



lS\, (5) 



with expansion 

e = -sign(y,)^, (6) 

and acceleration 

where i = 1,2,3, x' = {x,y,z}. Note that (||) differs from the expression usually given for 
the expansion (in Krasinski 1997, equations (4.1.4) and (4.9.6), for example) by the inclusion 
of the — sign(V|t) factor. Neglect of this factor is inconsistent since F enters the metric only 
quadratically, so the sign of the expansion cannot depend on the sign of F. The sign of 6 does 
depend on that of Vt, though: for the Fricdmann subcase of the Stephani-Barnes models, for 



example, V is related to the scale factor R{t) by V = 1/i?, and |i^V^t/V^| = 1 (see § |2.2|) , so 
that 9 — 3R/R = —SYt/V — ~3sign(l/f)/|i^|. This is important here because the constraint (0) 
contains the expansion. Note, too, that F is not a true degree of freedom parameterising distinct 
spacetimes, but rather represents a coordinate freedom, corresponding to different choices of the 
time coordinate. 

From (P), (^ and (1^), the condition (0) leads to the constraint: 

which is satisfied if and only if the function V has the form 

V{t, X, y, z) = T{t)X{x, y, z) + Y{x, y, z), (9) 

where T, X, and Y are arbitrary functions. This equation is the key additional constraint on the 
Stephani-Barnes solutions. 

It is worth noting that it follows from (||) that the acceleration scalar is constant along the fluid 
flow for every observer, i.e., ii^t = (where li^ = iiau"'), as can be seen by calculating -u^^ from (Q). 
In fact, it can be verified more generally that for any conformally stationary spacetime (i.e., a 
spacetime satisfying (m) and (0)), even with rotation, the acceleration scalar evolves according to 

where V denotes the spatially projected gradient (see Ellis 1998). It follows from equation (32) of 
Ellis (1998) that Vb9 — whenever the rotation vanishes (for a perfect fluid), which is the case for 
the Stephani-Barnes models. 

Note that if V has the form (0) we can immediately write the metric in a manifestly conformally 
static form 

v2 

ds^ = —J {-dT^ + X-^{dx^ + dy^ + dz^)} , 

where dr = T^tFdt, which shows that these models will indeed be conformally stationary, as 
required by Theorem P. We now discuss each subcase in turn. 



2.1. The Barnes Solutions. 

The Barnes solutions (Barnes 1973) all have spherical, plane or hyperbolic symmetry 
(i.e., they possess three-dimensional isometry groups acting on two-dimensional orbits, cf. i |2.2| ). 
The restrictions on the metric function V depend on which of these symmetries the spacetime 
possesses (see Krasihski 1997). For the solutions with spherical symmetry (the Kustaanheimo-Qvist 
solutions) or planar symmetry we introduce a new independent variable u{x,y,z) defined by 
u = r'^ = x'^ + y'^ + z'^ in the spherical case and u = z in the planar case. Then V is defined 
hy V = V{t, u), subject to the condition: 

= /(.)F^ (10) 

where / is an arbitrary function. Since V = V{t, u) we know from (S) that in order to admit an 
isotropic radiation field for all observers V must have the form 

V{t,u)=T{t)X{u) + Y{u). (11) 

Equation ( |lO|) then imposes a constraint on the functions X and Y, which will be outlined below. 

For the solutions with hyperbolic symmetry the constraint on V is very similar. This time we 
introduce the variable u = x/y. V can then be written 

V{t,x,y)^vW{t,u), (12) 

where W satisfies 

-Q^-fi^)W\ (13) 

with / once again a free function. The condition (^ now gives 

V{t, X, y) = T{t)X{x, y) + Y{x, y) = yW{t, u). 
Dividing by y and redefining X and Y in an obvious fashion we obtain 

W{t, u) = T{t)X{u) + Y{u), (14) 

in which X and Y are again constrained by dlq). 

For all three symmetries of the Barnes solutions, then, the constraints on the metric function V 
essentially reduce to the second-order differential equations (HQ) or (13). Imposing the condition (O) 
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) or (PJ) into (hm or (h 



introduces the additional constraint ( |ll| ) or (14). Substitution of (1 

respectively and differentiating twice with respect to time, dividing by Tt each time (and recognising 

that Tt 7^ 0, X ^ 0, so that F* ^ in (||)), leads directly to the condition 

f{u) = 0. 

Barnes solutions with f{u) = are conformally flat (Krasiiiski 1997, p. 142) and are therefore 
actually a subcase of the Stephani models. That is, proper Barnes spacetimes can be ruled out: 
they do not admit an isotropic radiation field. It only remains to apply the condition (H) to the 
Stephani models, which we do in the next section. 



2.2. The Conformally Flat Solutions 

The conformally flat sub-case is the entire class of conformally flat, expanding, perfect fluid 
solutions, and is the Stephani solution (Stephani 1967a, b). The function V is most often written in 
the form (see Krasihski 1997; Barrett and Clarkson 1999a, b): 

Vit,x,y,z) = ;^ (1 + i^Wlx - xo(i)n (15) 

In general the five functions of time in V are free. For our purposes, however, it turns out to be 
more convenient to use 

V = a{t)+b{ty ~2c{t)-r, (16) 

as in Barnes (1998), again with five free functions (we adopt three-dimensional vector notation, 
so that c = (ci,C2,C3), for example). In fact, this form is slightly more general than (^5|) - see 
Barnes (1998). 

We must be able to write V in the form (^ for the spacetime to admit an isotropic radiation 
field for all fundamental observers. From ( p^ or ( p^ ) it is clear that the functions X{x^), and Y{x^) 
in (0) can be at most quadratic in the x*. Writing X, and Y as quadratics in (O) and equating 
in dig ) all powers of x\ we obtain the following constraint equations: 

a{t) = aiT{t) + a2, 

b{t) - feir(t) + 62, (17) 

C(i) =CiT(t)+C2, 

where T{t) is a free function of time and the ai.2, 61,2 and Ci^2 are ten independent constants. Not 
all of ai, &i and Ci can be zero (in order that V^t 7^ in (^). 

Equations (p7|), along with (|4|) and (|lj), provide the complete set of irrotational perfect fluid 
spacetimes admitting an isotropic radiation field. Not all of the possible choices of parameters give 
rise to distinct spacetimes, though, and we outline in appendix ^ how coordinate transformations 
may be used to eliminate many of the parameters in ( |l7| ) , and determine when the models can be 
reduced to manifestly spherically symmetric form (c = 0). Given the forms for a, 6, and c in (O) 
we can equate ( |l5| ) with ( [l^ ) to find the corresponding constraints on i?, fc, and Xq in (15). This 
we do in appendix 0. 

At this point we can say that perfect fluid spacetimes admitting an isotropic radiation 
field for all fundamental observers are FLRW if and only if the acceleration of the fundamental 
observers is zero. This follows because the Stephani models with zero acceleration are FLRW (see 
Krasinski 1997, although it can easily be seen from (0): if m = then V ~ T{t)X{x'^) for some 
functions T and X, showing that a, b and c depend on the single free function T and so must be 
FLRW by the results of the next section). Thus we have proved: 

Theorem 3 The irrotational perfect fluid spacetimes admitting an isotropic radiation field for the 
fundamental observers are Stephani models with the free functions restricted by ^ITj) (or ^A^). 
These spacetimes are FLRW if and only if the acceleration of the fundamental observers is zero. 

It is worth noting that all of the models admitting an isotropic radiation field are manifestly 
conformal to (part of) the Einstein static spacetime (c/. Barrett and Clarkson 1999a) once they 



have been transformed so that c = cz as outlined in appendix ^. From ( [iq ) and (AS) we obtain 

V,t = a^Tt (l + — r^ 

since c^f — when c = cz (as noted in appendix K^ we can assume ai 7^ 0). Changing the time 
coordinate via dt ^^ aiT^tF dt, the metric (|j) becomes 

, (1 + jAr^f f 9 1 9 9 9^ 

where A = Abi/ai. The factor in braces is the Einstein static metric. Barrett and Clarkson (1999a) 
used this conformal relationship to simplify the study of the observational characteristics of a 
subset of these models. 



2.3. Symmetry and Thermodynamic Schemes. 

Having obtained the conditions (y_7|) for a Stephani model to admit an isotropic radiation field 
it is possible to say immediately that all such spacetimes possess (at least) a three-dimensional 
symmetry group acting on two-dimensional orbits (just as for the general Barnes models). This 
follows from the work of Barnes (1998), who showed that the dimension of the isometry group of 
any Stephani spacetime is determined by the dimension d of the linear space spanned by the five 
free functions a, b, and c: 

1. if d = 4 or 5 (i.e., at least four of the free functions are linearly independent), then the 
spacetime has no Killing vectors; 

2. if d = 3 there is a one-dimensional isometry group; 

3. if d = 2 there is a three-dimensional isometry group acting on two-dimensional orbits; 

4. if d = 1 there are six Killing vectors and the spacetime is Robertson- Walker. 



It is clear from (17) that a, &, and c depend on (at most) only two functions of time: /i(i) = T{t) 
and the constant function f2{t) = 1 (since T^f ^ we must have Tt 7^ 0, so that these are necessarily 
linearly independent). Thus, d — 2 and the solutions have three-dimensional isometry groups as 
claimed. If 02 = ^2 = C2 = in (O) then d = 1 and the spacetime is FLRW (see appendix ^. 

It follows further from this and the work of Bona and Coll (1988) (see also Krasiiiski, 
Quevedo and Sussman 1997) that all Stephani models that admit an isotropic radiation field for all 
fundamental observers also admit a strict thermodynamic scheme (that is, entropy and temperature 
functions can be found that depend on the energy density and pressure and satisfy the first law 
of thermodynamics as embodied by the Gibbs equation). The converse is not true, however, since 
there are Stephani spacetimes with d = 2 (which must admit a thermodynamic scheme) that cannot 
have an isotropic radiation field (these are models for which the second independent function /2(i) 
is not restricted to be 1). So, we have the following corollary to Theorem ||: 

Corollary 1 An irrotational perfect fluid spacetime that admits an isotropic radiation field has 
spherical, planar or hyperbolic symmetry and admits a strict thermodynamic scheme. 
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On the subject of thermodynamics, let us mention for completeness that the thermodynamic 
scheme occurring most often in the literature is that of a barotropic equation of state. It is known 
that the only Stephani models with a barotropic EOS are precisely the FLRW models (Bona 
and Coll 1988; Krasihski 1997). Thus, the only spacetimes with a barotropic EOS admitting an 
isotropic radiation field are FLRW models. This also follows (when fi + p ^ 0) from a theorem of 
Coley (1991). 



3. QCDM Models. 

The type la supernova results of Riess et al. (1998), Schmidt et al. (1998) and Perlmutter 
et al. (1999), which suggest that the expansion of the universe is accelerating, have lead to an 
increased interest in cosmological models in which a significant contribution to the energy density 
comes from either a cosmological constant or a scalar field (quintessence component), which is 
capable of driving the expansion (Peebles and Ratra 1988; Ratra and Peebles 1988; Caldwell 
et al. 1998; Zlatev et al. 1999). In QCDM models the matter is an admixture of non- interacting 
cold dark matter (CDM), i.e., dust, and a scalar field. The fundamental observers (galaxies) are 
implicitly identified with the geodesic congruence of the CDM. Note that there is no reason a priori 
why the scalar field gradient (which defines a natural 'velocity' field) should be aligned with the 
CDM congruence (although it will turn out that they are aligned if the fundamental observers see 
an isotropic radiation field). 

It is interesting to ask whether the EGS theorem can be extended to this case. We demonstrate 
that it can by proving the following theorem: 

Theorem 4 Any solution to Einstein's equations in which the matter consists of non-interacting 
radiation, expanding dust (CDM), and a scalar field (or cosmological constant), and for which the 
dust sees an isotropic radiation field, must either be an FLRW model, or have the gradient of the 
scalar field orthogonal to the dust congruence. 

(Note that the latter possibility means that gradient of the scalar field is spacelike, and is usually 
rejected as unphysical - although see below.) 

Proof: 

We may divide this proof into two parts: first we demonstrate from Einstein's equations in the 
1+3 formalism that any energy flux component with respect to the CDM frame must be zero if 
the CDM observers see isotropic radiation, then we show that the contribution to the energy flux 
(with respect to the CDM frame) from the scalar field is zero if and only if the gradient of the 
scalar field is parallel (or orthogonal) to the CDM velocity w" ^^ , so we deduce that the velocity 
fields are parallel (or orthogonal). The case where the field gradient is orthogonal to w"^^^ is 
probably unphysical, and will be rejected. Thus, the mixture of radiation, CDM and scalar field 
can be written as a single perfect fluid with geodesic fundamental congruence u° = u't^^^ , and it 
follows from the results of §0 that the model is necessarily FLRW. Throughout this section equation 
numbers of the form (E32) will refer to the lectures of Ellis (1998). 

Since the radiation is isotropic for the dust observers w" the energy-momentum tensor for 
radiation may be written in the perfect fluid form (0) with p = ^ji and the total energy-momentum 
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tensor is: 

Tab = f^UaUb + ^Hhab + pUaU^ + (f) ,a(l> ,b ~ 9ab {^<l>,c<t>''' + *(0)) , (18) 

-^ ^ ' "^v^ ^ « ' 

Radiation CDM Scalar Field 

where $(0) is the scalar field potential (often assumed to be zero, in which case the scalar field can 
be interpreted as a stiff perfect fluid). Note that the cosmological constant case can be included by 
setting (j) = A =constant, <&(</>) = 0. 

1. The fundamental congruence u°- is geodesic (u" — 0) because the CDM component does not 
interact with the other matter. This, in fact, implies that the rotation of u" must also vanish: 
from the momentum conservation equation for the radiation, we can write 

iia = --Valn^ = 
(where Va denotes the spatially projected gradient), so that (using (E27)) 

= V[aUb] = TVffcVa] In^ = -i^ab- = i;^bad , 

4 ^ ' 4 /i 3 

and we see that ujab = when 6^0. 
When iia and ujab are zero, (g) becomes 

V[a(6'Wb]) =UlbVa]d = 0. 

This implies (since Va^ — ^ ad — Sua) that 

VaO = 0. (19) 

(i.e., the expansion is homogeneous). From the constraint equation relating the divergence of 
the shear to other kinematical quantities (E32) we see that any energy flux component with 
respect to the CDM velocity field must vanish: 

qa = \^ae - 0. (20) 

This is the key step in the proof. 

2. Decomposing (n8) with respect to u° we find that the relative energy flux component is 

G = qa = -h^u'n, = -4>Va4>- (21) 

So (7a = if (/) = (the scalar field gradient is orthogonal to u°, and therefore spacelike), or 
if Va0 = (the scalar field gradient is parallel to u°). We take the latter case to be most 
important since the gradient of a scalar field is usually assumed to be timelike. 



Since Va0 = V a4> — Uacf) — ~Ua(j) it is possible to write (18) in the form of a single perfect 
fluid (0) with geodesic, shear-free, rotation-free velocity field; it is thus an FLRW model by §0. 

D 

It is easy to see from the above proof that the fact that the fundamental observers correspond 
to dust-like matter was not used, only that they followed geodesies. Hence the above result 
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applies for more general perfect fluids in place of the CDM component, as long as the fundamental 
congruence is geodesic. 

The idea of a spacelike scalar field gradient seems physically unappealing. However, such a 
field can (depending on the potential $) satisfy the weak, strong and dominant energy conditions. 
The strong energy condition will be satisfied if and only if $(</)) < everywhere, whereas as the 
weak and dominant energy conditions will be satisfied if <&(0) > 0, although not only so. Thus, a 
massless scalar field (stiff perfect fluid) with spacelike gradient satisfies all energy conditions. It 
should be borne in mind, though, that scalar fields arising in cosmological contexts often fail to 
satisfy the energy conditions. This case may deserve further consideration. As can easily be seen, 
the scalar field component gives rise to anisotropic stresses in the energy-momentum tensor, so 
such spacetimes are not FLRW. 



4. Conclusions. 

We have proved that the irrotational perfect fiuid spacetimes admitting an isotropic radiation 
field are Stephani models restricted by (|l^) (see also equations ([A^)), and are FLRW if and only if 
the acceleration u° of the fundamental congruence is zero (Theorem |^) . It follows from the fact that 
the constraints (O) depend on only two independent functions of time that all of the acceptable 
models possess three-dimensional symmetry groups acting on two-dimensional orbits (i.e., have 
spherical, planar, or hyperbolic symmetry) and therefore possess a thermodynamic interpretation. 
We have also shown that spacetimes containing a mixture of radiation, dust and scalar field (QCDM 
models) for which the dust observers see the radiation as isotropic must always be homogeneous 
and isotropic (Theorem 0) unless the scalar field gradient is spacelike and orthogonal to the CDM 
congruence - a possibility we reject as unphysical. This result also relies on the geodesic nature of 
the fundamental congruence. 

Crucial, therefore, to the proof of homogeneity and the verification of the cosmological principle 
is the non-acceleration of the fundamental observers. Despite the intuitive appeal of cosmological 
models in which the fundamental observers are associated with a dust-like matter component, it is 
unacceptable to simply assume that we are geodesic observers, especially when such an assumption 
is, in principle, testable. Acceleration leaves a characteristic dipole signature in the redshifts of 
nearby galaxies that may be detectable using galaxy surveys. The physical principle underlying 
this effect is easy to see. For a set of uniformly accelerated observers ('galaxies') in Minkowski 
space each observer will see other galaxies redshifted or blueshifted in a dipole pattern, with the 
blueshifted galaxies lying in the direction of the acceleration, because during the light-travel time 
between galaxy and observer the observer's velocity has increased relative the velocity at emission, 
so that the galaxies the observer is travelling towards are blueshifted, and those it is travelling away 
from redshifted. It also follows from this that the magnitude of the dipole increases with distance, 
simply because the light-travel time from more distant galaxies is larger. In a cosmological context 
this acceleration dipole must be added to other terms contributing to the redshift of nearby objects, 
in particular the expansion. The method of Kristian and Sachs (1966) and MacCallum and Ellis 
(1970) gives, for any cosmological model with fundamental congruence u", the lowest-order term in 
the redshift z as a function of distance r and direction e° (a spacelike unit vector orthogonal to u". 
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denoting the direction of observation) : 




z =. Hor fl - 


Uae" (Jabe''e^\ 


Ho Hq J 



0(r2), (22) 

where Hq = ^Oq is Hubble's constant and the last term in brackets indicates the quadrupole 
introduced by the presence of shear. In ( p^ ) r can be any cosmological distance measure (area 
distance, for example) because for small r all such measures agree to first order. Note that just 
as the monopole (expansion) term increases linearly with distance according to the Hubble law, 
so does the acceleration dipole. This is important, because it allows the acceleration dipole to be 
distinguished from any dipole resulting from the peculiar velocity of our galaxy with respect to the 
cosmological average rest frame (usually identified with the CMB frame). Equation ( p2| ) applies in 
this rest frame, and any peculiar motion results in a doppler shift for each galaxy, which introduces 
an additional dipole component into the galaxy redshifts. This dipole is just a constant depending 
only on the peculiar velocity of our galaxy. A boost to the 'correct' rest frame can eliminate this 
constant component, but cannot remove the acceleration dipole because it is distance dependent. 
It is important to note in this context that the acceleration referred to here is not the same as the 
'acceleration dipole' resulting from the gravitational attraction by the Great Attractor overdensity, 
which is often calculated using galaxy surveys (see Schmoldt et al. 1999). 

Galaxy surveys are often used to measure a possible bulk flow of the local universe, that is, 
the difference, if any, between the rest frame of the local universe and the CMB frame, which 
in standard cosmological models should be the same (see Willick 1998). A simple extension of 
these techniques (Clarkson, Rauzy and Barrett, in preparation) permits the acceleration to be 
constrained by observations. However, preliminary results suggest that the constraints on u"" are 
quite weak: it appears not to be possible to conclude definitively that we are geodesic observers. 
The accuracy of u°' determinations is limited both by uncertainties in the distance estimates to 
galaxies as well as the peculiar velocities of galaxies. 

Even if the acceleration was measured to be zero, it is still necessary to show that there are 
no anisotropic stresses (Ferrando et al. 1992) before the cosmological principle can be verified. It 
follows from equation (31) of Ellis (1998) that this is equivalent to determining that the electric 
part of the Weyl tensor is zero. It is not clear how this may be achieved using observations. 

Of course, the Copernican principle (which is a vital element of EGS-type theorems, allowing 
the high isotropy of the CMB here to be assumed for other points in the universe) remains a purely 
philosophical assumption. It has been suggested by a number of authors (see Goodman 1995, for 
example) that the Sunyaev-Zeldovich effect might be used to place constraints on the anisotropy of 
the CMB at distant positions, but it is not obvious that such observations will provide a definitive 
verification of the Copernican principle. Nevertheless, the arguments in favour of the Copernican 
principle are quite powerful, and it is a much weaker assumption than the cosmological principle. 
Note that if the acceleration is measured to be zero here, the Copernican principle must also be 
applied to give geodesic observers everywhere for the results of this paper (and the other ECS 
papers) to hold. 

Finally, one might expect that the 'almost' version of Theorem ^ would lead to spacetimes 



that are almost the Stephani models of (17). However, when the assumption of geodesic observers 
is relaxed it is no longer possible to constrain the rotation to be small, and the class of perfect fluid 
spacetimes with an almost isotropic CMB may include examples with distinctly non-zero rotation, 
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unless other constraints are brought to bear. ActuaUy, it is possible that all spacetimes admitting 
an IRF are irrotational (see Coley 1991, §2.2), although this has not been shown definitively. It 
would be interesting to determine the class of all perfect fluid spacetimes admitting an isotropic 
radiation field, and if it turns out that they are indeed all irrotational then the Stephani spacetimes 
defined by (p^), (^ and (jl^) are indeed the complete set. 
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A. Coordinate transformations and the Stephani Spacetimes. 

On the face of it the Stephani models admitting an isotropic radiation field in (|l^) depend on 
one free function and ten free parameters. However, it is possible to use coordinate transformations 
on the spacetime to eliminate many of these parameters, resulting in a considerable simplification. 
As is shown in Barnes (1998), conformal transformations of the coordinates on the hypersurfaces 
of constant time preserve the form of the metric but change the free functions a, 6, and c. These 
transformations can be thought of as acting on the five-dimensional space spanned by a, b and c 
and constitute the Lorentz group in five dimensions, 50(4, 1): they leave —ab+ |cp invariant (a 
and b are 'null coordinates'). It will be convenient here to let a = a + (3 and b — a — j3 (so that the 
transformations preserve —a^ -1-/3^ + |cp), and to adopt five-vector notation: q^^ — (a,/3,c). We will 
use the terms 'rotation' and 'boost' to refer to the transformations on g, and will call q timelike, 
spacelikc or null if — a^ +/3^ + |c|^ is negative, positive or zero, as usual. Then it is easy to visualise 
the transformations on the free functions by imagining the 'mass hyperboliods' of representations 
of the Lorentz group in the usual way: a timelike vector can always be boosted so that it has the 
form (a, 0,0), whereas a spacelike vector can be boosted and rotated into (0,/3, 0), for example. 

In addition to the Lorentz transformations we also have the freedom to change basis in the 
function space spanned by the free functions. For the spacetimes of interest here, described by (|l7|), 
it is desireable to preserve /2(i) = 1, so that the basis change is 

T^jT + S. (Al) 

In five-vector notation the equations ( p7| ) become 

q'^{t) = {a,P,c)^q'^T{t) + q^ (A2) 

(with qi and (72 constant vectors), and the goal is to reduce as many of the components of g^ 
and ^2 to zero as possible using Lorentz transformations in the five-dimensional space containing 



qi and (72, and the basis change (Al). Note that the FLRW (d = 1) subcase of ( |17| ) is characterised 
by the linear dependence of gi and 52 ■ 

It is easy to see that c (which breaks the spherically symmetry of the metric (^)) may always 
be reduced to the form c = cz (z = (0,0, 1)), with c a constant: perform a spatial 4-rotation to 
reduce g^ to Qi = (ai,/3i,0), then a spatial rotation amongst the c-components (which obviously 
leaves gi unaffected) to give g2 = (a2,/32,cz). 

It is possible in general to make further simplifications, but precisely how gi and g2 are 
simplified depends on whether they are spacelike, timelike or null. For example, if either gi or g2 



is timelike (or may be made timelike by a transformation (Al)) it is possible to reduce the model 
to manifestly spherically symmetric form (c = 0): boost so that the timelike vector, say gi, 
becomes gi = (ai, 0, 0) and rotate spatially so that the c-components of the other vector are also 



zero, g2 — (a2,/32,0) (we could then use (Al) to eliminate more of these constants). 



To summarise, we have demonstrated that it is always possible to reduce the Stephani models 
of (Irn) to the form 



(A3) 



a(i) = 


= aiT(t) + 02 


m-- 


= biTit) + 62, 


c(i) = 


= cz. 
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(where we have transformed back from a and (3 to a and b) , and when either of the qi or (72 is (or 
may be made) timehke we can set c = 0. 



Finally, note that we may always assume that ai 7^ in (|A3| ), because if ai = then 61 7^ 0, 
otherwise V^t = 0, and it is possible to perform a coordinate inversion x i~^ x/r^ that interchanges 



a and b). This does not exhaust the possibilities for simplification: we could, for example, use (Al) 
to set a2 = 0. 



B. The Constraints on R, k and xq. 

To find the constraints on R, k and xq in (nsl) corresponding to ( p7| ) first equate powers of : 
in (15) and (p^) to obtain 



(Bl) 
(B2) 
(B3) 

(B4) 

(B5) 

(B6) 

which are valid whenever ab — jcp 7^ (otherwise V cannot be written in the form (pa)). 

To impose the constraints ( |l7J ) perform the transformations of appendix |A| so that c = cz as 











a{t) = 


R^ 4R 


Xo 










m = 


k 
4i?' 












c(t)- 


k 
4i?"°- 




Solving these 


equations 


for R, 


k and 


xo gives 
R{t) 


b 






ab — 


C|2 










\Ht) 


62 






ab — 


C|2 












c 












Xo 


" 6' 





in A3) and ai =/= 0. Then |c| = c , and 6 and a are related by 



b — — a - 
ai 



62- 






7a + S, 



(B7) 



where 7 = bi/ai and ^ = 62 — 6102/01 are constants. Using this in (pj) and (B5) leads, after some 
rearrangement, to a quadratic relationship between k and i?: 

2 



(7 + <5i?) 



7c^i?^ = 0. 



In addition to this constraint relating k and R we can trivially rewrite (B3) or (B6) as 



4i? 
xo = -cz. 



(B8) 



(B9) 



Equations (pq) and (B£) are the constraint equations on R, k and xo corresponding to the 



equations (17), or rather (A3). If desired (B8) can be solved to obtain 

k _ 1 
4^2 



j + dR± V(7 + '5i?)2+47c2i?2 



- 18- 
From (B8) it is clear that the spherically symmetric Stephani spacetimes adraitting an isotropic 

-ij + SR)) =0, 



radiation field, for which c — 0, satisfy 

k fk 



4 V4 
which has the solutions fc = and R{t) free (flat Friedmann), or k linearly related to R, 

\k{t)=-i + 5R{t), 

with R{t) again free (when 5 = Q these become Friedmann models with curvature k — 47). The 
Stephani models studied by Barrett and Clarkson (1999a, b) are members of this class (with R 
depending quadratically on t), which explains the isotropy of the microwave background found for 
those spacetimes. 



